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ABSTRACT 

The Pan-STARRS 1 (PSl) survey of M31 (PAndromeda) is designed to 
identify gravitational microlensing events, caused by bulge and disk stars (self- 
lensing) and by compact matter in the halos of M31 and the Milky Way (halo 
lensing, or lensing by MACHOs). With the 7 deg^ FOV of PSl, the entire disk 
of M31 can be imaged with one single pointing. Our aim is to monitor M31 
with this wide FOV with daily sampling (20 mins/day). In the 2010 season we 
acquired in total 91 nights toward M31, with 90 nights in the rpi and 66 nights 
in the i-pi. The total integration time in rpi and ipi are 70740s and 36180s, 
respectively. As a preliminary analysis, we study a 40' x 40' sub-field in the 
central region of M31, a 20' x 20' sub-field in the disk of M31 and a 20' x 20' 
sub-field for the investigation of astrometric precision. We demonstrate that the 
PSF is good enough to detect microlensing events. We present light curves for 6 
candidate microlensing events. This is a competitive rate compared to previous 
M31 microlensing surveys. We finally also present one example light curve for 
Cepheids, novae and eclipsing binaries in these sub-fields. 
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Introduction 



Since the seminal paper by iPaczynskil (119861 ) . microlensing has been used as an el- 



egant way to probe the massive compact halo object (MACHO) content in our Milky 
Way. Various experiments have been conducted (e.g. MACHO, EROS, OGLE, MOA) 
targeting the Magellanic Clouds to constrain the MACHO mass fraction in the Milky Way 
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On the other hand, ICrottsI (119921 ) proposed to target M31 to search for MACHOs as 
dark matter candidates in spiral galaxies. This is because we can have multiple lines of 
sight toward M31 to better discriminate self-lensing from halo lensing. For example, the 
far-side of M31 experiences a higher column density of MACHOs than the near-side, thus a 
spatial asymmetry of the lensing events should show up if a significant MACHO component 
exists in the halo of M31. A description of the current statu s of M31 micro l ensing searches 
and the theory of microlensing in crowded fields is given in ICalchi Novatil (120101 ) and ref- 
erences therein. To summari ze this review, pre vious M31 projects (e.g. POINT-AGAPE 
jAuriere et a/.ll200lh . MEGA JCrotts et a/.ll200lh ) have monitored ^ 30' x 30' large fields of 
M31 on the near an d far side of the disk while avoiding the central bulge. The WeCAPP 
(IRiffeser et a/.ll200ll ) team took another approach and monitored a 17.2' x 17.2' field centered 
on the M31 center. They have detected ~ 10 events in their 11-year campaign (Riffeser et 
al. in prep.) 

The result of t he POINT- AGAPE team suggests the existence of MACHOs to explain 
some of the events ( Calchi Novati et al. 20051). whil e the MEGA team claims that the result 



of their 14 r nicrolensing events (Ide Jong et a/.ll2006l ). though 2 of them are more likely to be 
supernovae (ICseresnjes et a/.ll2005l ). is reconcilable with self-lensing s cenario. Few possible 



issues on the MEGA results have been addressed by the analysis of llngrosso et al\ (12006 
20071 ) . One aspect is that the dust in the M31 disk jMontalto et a/.ll2009r yields an asymmet 



ric detection efficiency of true variables and microlensing events and thus has to be included 
in a quantitative modeling. The other aspect is that the total number of events found is still 
not large enough to allow for small errors in statistical analyses on the MACHO fraction 
of M31. The investigation of the asymmetric signal, the spatial distribution and the total 
number of events requires large fields to be observed. Another point is the time sampling: 
Because the full-width-half- maximum (FWHM) timescale of the mic rolensing events in M31 



are estimated to be shorter than a few days (IRiffeser et adl2006l ). one requires a higher 
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focused on the 17.2' x 17.2' FOV in the bulge of M31 with daily observations. 



Ongoing campaigns attempt to monitor the inner region of M31 with higher cadence. 
For example, the ANGSTROM (Andromeda Galaxy Stellar Robotic Microlensing) collabo- 
ration employs a world-wide newtork of 2-me ter telescopes to al low for 24-hour survillance 
of the inner bulge (4.6' x 4.6' FOV) of M31 JXerins et a/.ll2006h . The PLAN (Pixel Lens- 
ing Andromeda) collaboration aims for consecutive-nigh ts campaigns for two field s (with 
13' X 12.6' FOV each) around the inner region of M31 (ICalchi Novati et a/.l 120091), which 
lead to the first finite- source microlensing event, 0AB-N2, toward M31( Calchi Novati et al. 



2010). The analysis on this even t, as well as the bright event GLl /PA-S3( Riffeser et al. 20081 ) 



and PA-Nl (jAuriere et a/.ll200ll ). favors a MACHO senario relative to self-lensing. However, 
a comprehensive study of the MACHO content in the halo of M31 and constraints on the 
self-lensing rate requires a large FOV survey with daily sampling. This is only possible with 
the advent of PAndromeda. 

PAndromeda is one of the 12 key projects within the scope of Pan-STARRS 1 (see 
Kaiser et a/.ll2010l : iHodapp et a/.ll2004l : iTonry fc Onakall2009l . for a detailed description of the 
PSl system, optical design and the imager). PAndromeda is designed to identify microlensing 
events toward M31 with high cadence observations (0.5 h per night for a time span of 5 
months per year). With the large FOV of Pan-STARRS 1, it is possible to monitor the 
entire disk of M31 with one single exposure (see Fig. [1]). This enables us to compare the 
self-lensing event rate (cau sed by bulge and d isk stars within M31) to the prediction from 
the theoretical calculation ( Riffeser et al. 20061 ) with a certain detection efficiency (which is 
given from the set-up of PAndromeda and simulations). The merit of PAndromeda is the 
ability to perform a differential measurement across the entire M31 disk. With similar PSF 
and exact cadence for each image and filter, the detection efficiency mostly varies due to 
M31 background light. The impacts of PSF variations and cadence are rather small. A 
discrepancy between the observed and theoretical self-lensing event rate can indicate that 
the assumptions in the M31 model under consideration (e.g. stellar population) might need 
to be improved flSaglia et a/.ll2010l ). 



The microlensing event rate will help us to constrain the mass-fraction of compact 
objects in M31 and the Milky Way halo, as well as the mass function at the low mass end of 
the M31 bulge. PAndromeda further aims to shed light on the stellar population properties 
of M31 based on the color profiles, surface-brightness (SFB) fiuctuations, resolved stars and 
variables. It will thus improve our understanding of the mix of stellar ages and metallicities 
in the halo, bulge, stellar streams and dwarfs of M31. This information is required for an 
accurate interpretation of the microlensing events. For long timescale events, PAndromeda 
has the potential to detect sub-stellar objects by means of excursions from the standard 



light curve due to co mpanion objects (similar to the candidate discussed by lAn et a/.ll2004b 



Ingrosso et a/.ll2009l ). In addition to microle nsing, the high-cadence of PAndro meda data-set 
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This paper is organized as follows. In section [2] we present the survey parameters 
and observations done in 2010. We give a detailed description of our data reduction, data 
analysis and demonstrate the data quality of the 2010 observation season in section [3J The 
variables detected in the PAndromeda data-set are presented in section |H An overview of the 
microlensing events is shown in section [5l followed by an outlook in section |6l We conclude 
the paper in section [3 



2. PAndromeda Survey Parameters and Observations in 2010 

The PAndromeda project is one of the 12 key projects of PSl. It monitors M31, the 
Andromeda galaxy. The PAndromeda survey is also one of the Pan-STARRS medium deep 
fields (MDFs) and is called MDO alternatively. PSl observations are carried out with the 
1.8m Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) located at 
Haleakala in Hawaii. The camera used is currently the largest one in the world. 

It consists of 60 detectors in an 8x8 array with empty corners. Each detector is seg- 
mented in an 8x8 array of 590x598 pixel cells with gaps between cells of 19 and 13 pixel 
widths in the column and row dimensions, respectively. In addition, between each detector 
there are gaps of approximately 300 pixels in both directions. A single detector therefore 
consists of a 4872x4824 pixel array. Since pixels have a size of 10 /im which corresponds to 
0.258" in the focal plane the FOV of each detector is 20.95' x 20.74' and the total FOV of 
the 60 detector array is about 7 degrees^. The layout of the detector surface is illustrated 
below in Fig[T] 

PAndromeda started the first observation season in 2010 and monitored M31 from 
23/07/2010 to 27/12/2010. The M31 observations are done in rpi and ipi-band. We use two 
filters to confirm microlensing events from their achromaticity. We acquired 90 nights in rpi 
and 66 nights in ipi. The total number of images and integration times obtained in the first 
observing season of PAndromeda are listed in Table [U 

The 0.5 h integration time per night was divided into two observation blocks separated 
by 3 to 5 hours in order to trace microlensing events with timescales shorter than 1 day. The 
observation cadence is shown in Fig. [2j 
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Fig. 1. — The Giga Pixel Camera (GPC) and its 60 detectors labeled by the blue and red 
grids. The red numbering indicates the same detector corresponding to the blue numbering, 
but rotated with 90 degrees during the observation . Each detector is composed of 4872x4824 



pixel with 0.258" per pixel ( iTonry fc Onakall2009[ ). The field of view of GPC is ~ 7 degree 



^2 



The entire Andromeda galaxy, M32 and NGC 205 can be imaged with one pointing (as shown 
in the figure). The underlying M31 image is a single rpi-band frame, which was chosen to 
have relatively small masked area. The masked areas (in black) visible e.g. in the detector 
14 (marked in blue) are due to ill-functioning orthogonal transfer arrays (OTAs), video guide 
star data and areas on the focal plane with sub-standard imaging performance within this 
detector. These gaps can only be closed by large dithering or by rotating the camera. 
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Table 1: PAndromeda integration times in rpi- and zpi-bands for the 2010 season 



Filter 


Nights* 


Images 


Total integration time 


rpi 


90(74) 


1179 


70740s 


ipi 


66(56) 


603 


36180s 



*The number of nights with 2 visits are indicated 
in the parenthesis. 



From the data header information one can see that each of the two nightly observing 
blocks yielding 12 frames with 60 seconds of exposure in 2 filters needs 16.2 minutes of 
telescope time. The overhead of about 35 percent includes the 13 seconds of read-out time, 
15 seconds for filter change and about 30 seconds for focusing after a filter change. Due to 
the overhead for filter changes we requested to take only data in the rpi and ipi filters in 
2010. The rpi and ipi filters are chosen because the source of the microlensing events are 
expected to be red evolved stars and the contamination from intrinsic variables in M31 is 
less severe in the rpi filter. 

Note that the central part of the camera close to the bulge of M31 is slightly out of 
focus (see Fig. [1]). 

For each filter, the observations within one visit have a polygonal dithering pattern with 
a radius of 9". It is either a hexagon (rpi-band) or a square (ipi-band) with the pointings 
located at the edges and the center. 



3. Data processing 



All images obtaine d for the PAndr omeda project are processed through the Image Pro- 
cessing Pipeline (IPP) f lMagnierll2006[ ). The IPP runs the images through a succession of 
stages, including bias and dark correction , masking and artifact removal, flat-fielding, as- 
trometric calibration ( iMagnier et a/.l 120081 1 and a fiux-conserving warping (the images are 
resampled to a common pixelscale) to a sky-based image plane, so-called skycells. In the 
following we use the term "image" for the warped skycells. 

The warped images have a dimension of 6000x6000 pixels and a plate scale of 0.2"/pixel 
which is smaller than the natural plate scale of 0.258"/pixel of the PSl Giga Pixel Camera. 
Each image has an attached variance map and a 16-bit mask image. For the PAndromeda 
project all images provided by the IPP are further processed using the PAn dromeda data re- 
duction pipeline. This pipelin e is based on the difference imaging technique (jAlard fc Lupton 
19981 : iGossl &: Riffeserl 120021 ) which provides high precision photometric measurements in 
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Fig. 2. — PAndromeda observation cadence. The histogram shows the number of visits (0, 
1 or 2 times per night) on M31. The black hne indicates the lunar phase (1 for full moon 
and for new moon). The red line shows the distance of the M31 center from the moon as 
a function of observing date. A pause in the observations is visible in December when we 
were reaching the limit of the 2% PSl-time. 



highly crowded fields. The PAndr omeda pipeline has b een integrated by us into the Astro- 
WISE data management system (IValentijn et a/.l 120071 ) which has been developed to serve 
surveys with large data volumes. For a detailed des cription of the Astro WISE integration of 
our software we refer to iKoppenhoefer et al\ (1201 ll ) . 



In the following we describe each processing step that is performed by the PAndromeda 
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pipeline starting with the warped images until the creation of the light curves and the 
detection of microlensing events. Fig. [3] gives an overview of the different processing steps. 
These steps are performed independently for each filter rpi and ipi and for each skycell. 

At the end of this section we analyze the astrometric and photometric stability of the 
light curves (see section \3M and section [3.10p . 



3.1. Interface between IPP and the PAndromeda Pipeline 

By default the IPP subtracts the background in each warped image using a bi-linear 
interpolation of a model with coarsely sampled super pixels. The procedure works well 
in case of pure sky background. In the PAndromeda survey, however, the images contain 
also astrophysical background, i.e. the brightness profile of M31 which is one of the largest 
objects on the sky. In the PAndromeda pipeline the brightness profile is actually used later 
on. Therefore we requested to restore the background by the IPP that originally had been 
subtracted. 

In the first data release the skycell layout used for the processing of the PAndromeda 
data was not optimized. In particular it does not make use of the fact that the detector 
boundaries are more or less the same with respect to M31 in each image. In some months 
from now a reprocessing of the PAndromeda data will be done by the IPP with skycells that 
are aligned with the outlines of the detectors. 

Although not yet optimal, the processing used for this work is still good enough for 
the detection of microlensing events. In total we obtained 70 skycells with acceptable filling 
factor. We require at least 50% of each image not to be empty, to save disk space and 
computation time by avoiding images with large empty areas. One example, skycell 077, is 
shown in Fig. |H Each skycell contains pixels of 4 different detectors in the current layout. 
In order to save disk space and to speed-up processing we apply a 2x2 pixel binning to all 
images. In this way we are able to produce the light curves of a significant fraction of the 
total FOV of the first season data within a few weeks of computing time. In the future (i.e. 
when all previously taken PSl images will be reprocessed with an improved Image Processing 
Pipeline) we will switch to the natural plate scale of 0.258"/pixel. 

In the first processing step of the PAndromeda pipeline all images are registered to the 
AstroWISE system, a procedure during which the IPP variance map and IPP mask images 
are merged into a single file - a variance map with zero values for masked pixels. The IPP 
variance maps have correlated noise. Due to remapping and reduced pixel size this correlated 
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Fig. 3. — Overview of the processing steps performed by the PAndromeda data reduction 
pipehne for each skycell. 
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noise is not accounted for in our pipeline, but we empirically correct for this by rescaling the 
error bars of the PSF light curves, (see section [3T1) . 

The PSF FWHM of each im age is determined using the median of all FWHM JMAGE 
values mea sured by SExtractor (Bertin fc Arnoutslll996l ). In addition we apply a filtering 
algorithm ( iGossl &: Riffeserl |2002| ) to detect and mask cosmic rays based on their narrow 
FWHM. 



3.2. Creation of a Reference Image 

The difference imaging technique requires the creation of a reference frame which is a 
stack of the images with the best image quality. In our case we select the 30 best seeing images 
for each skycell excluding images with a high masking fraction and/or a low transparency. 
We create a reference image in rpi-band and ipi-band for each skycell. The individual 
processing steps are as follows: 

In the first step we photometrically align each of the 30 input images Jj to the image with 
the very best seeing Ii. To account for zeropoint differences we multiply each image with 
a scaling factor /j. Furthermore we add a model Bi to account for differences in the sky 
background: 

Ii = f,xl, + B, . (1) 

For a proper photometric alignment of the single images, it is very difficult to calculate 
the scaling factor using a set of reference stars. In order to make use of all pixels in the 
field of view and to deal with the high crowding, each scaling factor /j is calculated with 
the difference imaging technique: using all unmasked pixels of the image we calculate an 
un-normalized convolution kernel that would match the PSF of image Ji to image Ii and 
would simultaneously scale image Ji to the flux level as image Jj. We do not use the kernel to 
really convolve the images but rather calculate its inverse norm which delivers us the scaling 
factor fi. In this way all scaled images I- have the same zero point but the PSF remains 
unchanged. Note that because the sky background in image Ii is still differing from the 
sky background in image /i at this stage, we need to account for a sky background change 
while doing the calculation of the kernel. This is done with including a 1st order polynomial 
background in the minimization procedure. 

The polynomial background is very robust when calculating the scaling factor fi, how- 
ever, for doing a proper adjustment of the sky background we do not use it and calculate a 
bicubic spline model Bi which has 16x16 nodes with a separation of 183 pixel. The spline 
values at each grid point are calculated in the following way: first we apply a 61x61 binning 
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to the scaled images fi x /j (excluding all masked pixels) and subtract the binned images 
from the binned version of Ji. The subtracted images have a dimension of 49x49 pixel and 
are then convolved with a 3x3 median kernel resulting in an array of 16x16. These values 
are taken as the spline values at the grid points which are used to calculate the sky back- 
ground correction Bi. After the photometric alignment all images have the same zeropoint 
and sky background. 



Recently, iKerins et al\ ( l2010l ) found that image subtraction in M31 works better when 
first performing a careful photometric alignment prior to matching the PSF. Although de- 
veloped comple tely independ e nt, ou r procedure is very similar. The difference of our method 



with respect to iKerins et al\ (120101 ) is that we account for PSF differences when doing the 



photometric alignment. 

The reference image is a weighted stack of all 30 photometrically aligned images where 
the weight, «j, for each image is calculated based on the background noise (median of a in 
each pixel) and the PSF FWHM, as measured by SExtractor during the image registration: 

a, = l/(a, X FWHM,)' . (2) 

Before stacking we replace masked pixels in each input image with the pixel values of the 
other image which has the most similar PSF. The similarity of the PSF is determined using 
a set of isolated stars. If a pixel is masked also in the most similar image, we replace with 
the second most similar image and so on. The procedure to replace masked pixels results in 
a very homogeneous PSF in the stacked image which is essential for the following processing 
steps. Note that this method only works if the set of input images already has similar PSFs. 
Fig. H] shows an example reference image of skycell 077 in rpi-band. The median PSF 
FWHM of all 30 input images is 0.861" with RMS = 0.012". In order to obtain the reference 
fluxes of all resolved stars in the reference image we subtract a bicubic spline model and 
perform PSF photometry on each source in the background subtracted image. The source 
positions are taken from a SExtractor table that has been created with the parameters 
DETECT_THRESH=3 and DETECT_MINAREA=5. We iteratively subtract neigh boring 



stars to remove blending effects (similar to the program DAOphot, IStetsod 119871 ). The 
eference fluxes will be added at a later stage to the fluxes measured in the difference images 
in order to create the light curves of resolved stars (see section [3^ . 



3.3. Creation of a stack for each visit 



Each visit of M31 consists of (up to) 7 individual exposures in rpi-band which are 
obtained in a dithered mode. For some visits we also obtained (up to) 5 individual exposures 



- 12 - 




Fig. 4. — Reference image of skycell 077 in rpi-band. The skycell contains pixels from 4 
adjacent detectors with the vertical and horizontal masked area being the gap between them. 
Note that the layout of the skycells has not been optimized for the PAndromeda survey. In 
the future we will use skycells that are aligned with the detector borders. Each skycell image 
has a size of 20' x 20'. 
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in ipi-band (see section [2]). The single exposures of each visit are combined to a visit stack 
for each filter in order to improve the S/N and to eliminate masked areas. The procedure 
to create visit stacks is identical to the procedure described in the previous subsection. Also 
here, masked pixels in one image are replaced by the pixel values of the image with the most 
similar PSF. Note that although the number of input images is lower than in case of the 
creation of the reference images the requirement of finding an image with similar PSF is 
satisfied since all images of one visit are taken subsequently after focusing the camera in the 
beginning of the visit. 



3.4. Photometric alignment 

The single visit stacks are photometrically aligned to the reference image using the same 
procedure we described in section |3^ After this step the visit stacks and the reference image 
have the same zeropoint and sky background and differ only in depth and PSF. Note that 
this zeropoint has an arbitrary value and only relative fluxes are calibrated. An absolute 
photometric calibration is done a posteriori as described in section 13.61 



3.5. PSF alignment and image subtraction 

For each single visit stack we subtract a PSF aligned reference image and obtain a 
difference image . The PSF alignment is done using the standard approach proposed by 



Alard fc LuptonI (119981 ) with an optimal convolution kernel modeled as the superposition of 
a set of Gaussian kernel base functions which are modulated by a two-dimensional polynomial 
function in the x- and y-direction. We use the standard kernel model with 4 base functions 
having standard deviations of 9, 6, 3 and 0.1 pixel and a polynomial order of 2, 4, 6 and 
respectively. The kernel is normalized (flux conserving) and has a size of 51x51 pixel. The 
total number of free parameters is 50. In order to account for PSF variations across the 
skycell we divide the image into 8x8 sub-fields and calculate a conv olution kernel for each 



sub- field separately. Our procedure is similar to iKerins et all ( l2010l ). The same procedure 



is also used for the WeCAPP data-set (Riffeser et al. in prep.). 

The difference images contain only noise at the positions of constant sources. Variable 
objects such as microlensing events or eclipsing binaries are visible as positive or negative 
PSF-shaped residuals. Fig. |5] shows an example difference image of the bulge region of M31 
in which several variable objects are clearly visible. Note that in the bulge of M31 the stellar 
density is so high that a large fraction of the pixels contains variable sources. 
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Fig. 5. — Zoom-ins to an example difference image of skycell 078 in ipi-band. Left panel: 
Positive (white) and negative (black) residuals are visible at the position of variable sources. 
Right panel: The center of M31. The white circle indicates a distance of 5' (1.12 kpc) from 
the center of M31 (marked with the white cross). The vertical and horizontal masked area 
are the gaps between the detectors. The size of this image is 20' x 20'. The white box 
outlines the size of the zoom-in presented in the left panel. The flux cut-level of these two 
images are the same. 



3.6. Creation of the light curves 

The light curves are created using PSF-photometry on the difference images. The PSF 
is constructed using a set of isolated stars in the convolved reference image. The pixels of 
each PSF reference star are subdivided into a sub-pixel grid, shifted and combined. Since 
the stars are all on different sub-pixel positions, we obtain a smooth profile for the combined 
PSF with sub-pixel precision. The PSF is then used to measure the flux in each difference 
image after binning to the actual sub-pixel position of each star. 

As shown in Fig. 12] we produce two sets of light curves. On the one hand we create a 
light curve of each pixel which we use for the detection of microlensing events. These light 
curves only contain differential fluxes as measured in the difference images. On the other 
hand we create light curves of resolved sources both in M31 and in the foreground. For 
those light curves we add the flux as measured in the reference image (see section 13. 2p to 
the differential fluxes and obtain absolute fluxes. 
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3.7. Flux calibration 

The flux calibration requires an existing cat alog with similar filter as PSl. Thus we 
use the catalog of Data Release 7 from SDSS-II (lAbazajian et a/.ll2009l ). The rpi and ipi 
are rather close to the SPSS ana logs thus the overall fiux calibration should be insensitive 
to color terms f lTonry et a/.ll2012[ ). The data used consist of the M31 footprints Run 3356, 
3366, 3367, 6426 and 7210 with a total of 6061956 sources. We calibrate our instrumental 
magnitude {'minstr,PSF) in the reference image to the SDSS magnitude {rn^^gg) with the fol- 
lowing formula: 



(3) 



^SDSS — '^instr,PSF + ZPpsF 

= -2.51og(F,e/[PSF,e/]) + ^PpsF 

where the i^re/pSF^e/] is the fiux of a source determined with the PSF of the reference image 
PSF^e/- Note that the PSF profile and the amplitude is arbitrary and thus different in each 
skycell. This is because the reference images are created out of 30 images having the best 
seeing, which do not necessarily have the same timestamp. The PSF profile is measured with 
selected stars on the reference image. ZPpsp takes into account the difference between the 
SDSS magnitude and 'minstr,PSF and also depends on the exact PSF-construction (how many 
and which stars are used, which kind of data enter the reference image). The flux calibration 
is performed on the reference images instead of the individual image, because the reference 
images have the highest S/N. 

Each difference image is calibrated to have the same flux unit as the reference image 
and therefore 



m,^,, = -2.5\0giFref[PSFref] + AF,,ff[PSFd^ff]) + ZPpsF ■ (4) 

holds. This is very useful to determine the magnitude of variables. In an extreme case, the 
magnitude of a transient event which has not been detected before (Fre/[PSF,.e/] = 0) can 
be easily determined from the flux in the difference image alone. 

The relative calibration of PSl data with SDSS data in the ipi and rpi bands does not 
take into account color terms, but still is fairly accurate. The rea son for this is that these 
color terms are small in the rpi and ipi fllters (see lTonry et a/.ll2012l ). At a later stage we will 
nevertheless calibrate our PAndromeda data using PSl calibration tools and data directly. In 
case the reader is interested in the mag nitude transformat ion to the Johnson/Kron- Cousins 
BVRI system, we refer to the work by llvezic et al\ (120071 ): 



R — r SDSS 

I — 'i'SDSS = 



-0.0107{rsDss - isDss) 
-0.0307(r5Z)55 - isDssf 



' + 



0.005(rsD5S - iSDSs? - 0.2689{rsDss - ^sdss) - 0.154 
0.1163(r5D55 - tsDssY - 0.33Al{rsDss - isdss) - 0.3584 

(5) 
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with the root-mean-square scatter for residuals evaluated for all their sample stars to be 15 
millimag in R — tsdss and 19 millimag in / — isnss- 

Our single frame data are remapped twice: there is a remapping to the sky tessellation 
where the pixelsize is decreased. Later on when the single frame prereduced data from IPP 
(Sect. 3) get ingested into our pipeline the images are rebinned to twice the pixelsize (Sect. 
3.1). The remapping causes a correlated noise depending also on the interpolation method 
used. Whereas our pipeline uses per pixel error propagation for all steps that are done after 
ingestion into our pipeline, it does not account for correlated noise that comes from the 
resampling before. We therefore empirically rescale the errors in the lightcurves directly, 
such that the formal PSF flux errors match the scatter of the PSF flux measurements. This 
rescaling factor is (as expected) the same for every skycell and every filter. 



3.8. Microlensing event detection 



The microlensing event detection is performed on the light curve of each pixel with 
successive criteria. The applied criteria for detection are as follows. For a preselection, we 
require at least three successive flux measurements in the light curve to be 3a above the 
base-line in eith er rpi or zpi-band. We then perform a microlensing fitting with the follow- 
ing formula (see iRiffeser et a/.l 120061 ) 



cff 



12(t-to) 



1 



+ aQ + ai{t - to) 



where F^g is the effective flux and for high magnifications is similar to the flux excess 
to is the time of flux maximum, tpwHM is the full-width half-maximum event timescale. 



(6) 



(7) 



The term ao takes into account the shift of the base-line flux by a constant in cases for 
which there is a variable spatially close to the microlensing event and data of these variable 
phases enter the reference image. The term oi is to include the gradient in the base-line flux 
in case the base-line flux is influenced by a long period variable close to the microlensing 
event. We require the data point closest to the time of the fitted maximum to have a good 
PSF. 



We then filter out light curves with a x = ^/Xdol ("where Xdof total divided by 

the degree of freedom) larger than 1.4 in rpi-band. We also require the S/N of the closest 
flux measurement to the time of the fitted maximum to be larger than 9 in both rpi and 
zpi-band. 
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The light curve of a microlensing event consists of many noisy data points that scatter 
around the basehne if one excludes the small time interval (of order 2 x tpwHM) where the 
microlensing event is measurable. The frequency of the deviations from the baseline depends 
on their amplitude and they occur at a rate given by the statistics of the PSF-flux errors. 
These deviations should be spread over the survey time and not be clustered in some time 
intervals. A correlated deviation from the baseline can be a hint of a nearby variable source 
or a time correlated artefact of data reduction. We would like to identify such correlations in 
the light curves and sort out microlensing candidates that have a good light curve fit but a 
signature for a time correlated deviation from the baseline not connected to the microlensing 
event. 

The Xdof value of a light curve fit only tells how well a model describes the data averaged over 
all data points, but it is not sensitive of "how" the contributions are spread over the data 
points, i.e. the survey time interval. We therefore add another criterion that quantifies the 
temporal correlations of the model mismatch. We define the model mismatch or deviation 
vector di = fi — f{ti) where fi is the i-th flux measurement at time ti and f{ti) is the model 
flux at time ti according to equation ([6]). The number of data points is n. With the further 
definition of Sj = sign {di}, so that Sj = 1 or Sj = —1 holds, we introduce 

n-l 

K = SiSi+i . (8) 
1=1 

En is a measure of how strongly the deviations of the data from the model ar e correlated 



in time. This equation is known also as the energy of a Lenz-Ising model ( llsind 119251 ) 
with nearest neighbor spin interactions only. Si is then the spin, and the spin-spin coupling 
constant usually called Jij is equal to minus one in our case. In n euroscience the Is ing model 



is often used to describe the energy of an ensemble of neurons flTang et a/.ll2008l ) and this 
also our motivation for choosing such an extra term. Equation ([H]) quantifies (similar from 
our perspective) the "energy" for the deviations of data points from a model as seen by 
neurons. 

For a random process, i.e. if the deviation vector d has no time correlation, the expec- 
tation value of equation is zero. The distribution of i?„ is gaussian with a width of ^/n. 
In order to make this scatter of in equation ([8]) independent of the number of data points 
(similarly as the Xdoi independent of the degrees of freedom) we normalize equation ([8]) to 

E = EJV^ . (9) 

We empirically require \E\ < 3 for the full light curve and l-E^'^'^'^l < 3.5 for the 20 fiux 
measurements closest to the time of the fiux maximum. We allow the time correlation of 
the deviation vector to be larger around the time of maximum brightness. This is purely 
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empirical. It can be justified by the need to allow for time coherent deviations from the 
points-mass-point-source standard model in the presence of extended sources and lenses, 
and disturbed deflection potentials. The criteria above (which are summarized below in 
Table [2]) are not fine-tuned to pick some particular PSl microlensing candidates but they 
where developed and tested for the WeCAPP project (Riffeser et al. in preparation). 

If the energy criteria l-El < 3 and < 3.5 are fulfilled, we select the light curves 

with tpwHM shorter than 20 days. This is to avoid the misinterpretation of long period 
variables as microlensing events due to the short time span of the first PAndromeda season. 
Contamination from the long period variables is less severe for campaigns with longer time 
span, e.g. the WeCAPP project (Riffeser et al. in preparation). All the criteria used are 
summarized in Table [2J All light curves that pass the microlensing event detection criteria 
are examined by eye. We exclude light curves for which the best-fitted maximum to and 
the time-interval plus/minus tpwHM around the maximum are not in the time-span of the 
PAndromeda data-set (event time selection). We also exclude the light curves with the 
signature of variability outside the event timescale. The latter happens, e.g., when there 
are variables or artefacts close to the microlensing candidate that contribute some flux to 
the examined PSF light curve. The number of event candidates that are discarded after 
the event time selection is of the same order as the number of our final events. Thus our 
detection method is mostly automatic. 

Table 2: Microlensing event detection criteria 
Criterion 

I Three successive 3a flux measurement in rpi or ipi 

II Good PSF around to 

III Xdof < 1-4 in rpi 

IV S/N > 9 for maximum flux measurement in rpi and ^pi 

V l-El < 3 for the full light curve 

VI l-^P'^^'^l < 3.5 for measurements around to 

VII tpwHM < 20 days 



3.9. Astrometric precision of the warped images 

In order to achieve high photometric accuracy with the image differencing technique 
it is very important to have a good alignment between the reference image and the single 
images - in other words to have a good relative astrometric calibration. In this section we 
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study the astrometric precision of the PAndromeda warped images and check the quahty of 
the astrometry achieved by the IPP by cahbrating against 2MASS. 

As a first test we compare the x- and y-positions of sources in the rpi- and zpi-band 
reference images of skycell 068 (see Fig. |6]). The RMS of the residuals is 175 mas (note that 
the distribution of residuals shows broad wings, and therefore the RMS value is dominated 
by large residuals. The standard deviation of a Gaussian fitted to the core of the distribution 
of residuals would be much smaller but this estimate strongly depends on the fitting radius). 
Further we check for systematic trends of the residuals as a function of position within the 
skycell (Fig. [7]). There is no trend visible proving a homogeneous astrometric calibration. 

The reference images have been constructed using the images with best image quality. 
In order to test if the internal astrometric calibration is good also for the rest of the data 
we compare the x- and y-positions of sources in two individual rpi-band visit stacks, one of 
them having a very low airmass (1.07) and the other one having a high airmass (1.54). We 
find a somewhat lower precision (RMS of 237 mas) but no trend with position within the 
skycell. 

We also tested the absolute astrometric precision comparing the positi ons measured in 



the r pi-band reference images of skycell 068 with a catalog from SDSS DR7 (lAbazajian et al. 



20091 ). Here we see a systematic trend only at the sub-pixel level in the RA coordinate as 
a function of RA as well as a constant systematic offset of the DEC coordinates (see Fig. 
IH]). We cannot say whether this points to a problem in the astrometric calibration done by 
the IPP or if the SDSS catalog has a calibration issue. In any case, for the quality of the 
light curves only the relative astrometric precision is important and since the offset is small 
compared to the pixel size we see no problems with identifying objects when cross-matching 
with external catalogs. 



3.10. Photometric stability and quality of the light curves 

To understand the data quality of PAndromeda, we first investigate two sub-fields, one 
in the bulge of M31 (skycell 077) and one in the disk (skycell 091). A histogram of the PSF 
distribution is shown in Fig. [91 The maximum peaks at 1". 

To quantify the precision of the photometry, we perform a test on sub-fields in the bulge, 
disk and the astrometric fields. We use SExtractor to extract the position of the resolved 
stars in the reference images of the sub-fields and create the light curves, i.e. flux difference 
of these stars relative to the mean flux as a function of time, at the position of these stars 
with our PAndromeda pipeline. Since most of the stars are not variable, the scatter of 
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Fig. 6. — Internal astro metric precision measured by comparing the positions of all sources 
in the rpi- and zpi-band reference images of skycell 068. 



the estimated flux difference (relative to the mean flux) is an estimate for the photometric 
precision as a function of the magnitude of the star. The ratio for this RMS-error relative 
to the stellar flux is plotted for each star in Fig. [TOl In the plot we present the 3a clipped 
RMS value. The theoretical S/N predictions, account for the contributions from sky, read 
out, object and scintillation noise, are shown in the black lines. We assume the surface 
brightness o f M31 in rpi to be 20.5 and 21.5 mag/arcsec^ at the location of skycell 077 and 
skycell 091 (ITempel et al \\2011\ ]. The RMS distributions from the data are shifted from the 
theoretical expectations. This is because masking is not taken into account. Besides, the 
resampling to the skycells and binning process also has influenced the quality of the light 
curves. Nevertheless, the result shows that the photometric precision is 1% for stars with 
rpi = 19 mag. 



4. Examples for variable light curves 

We also analyze the light curves of the resolved sources in the central 40' x 40' sub-field. 
We find about 60 Cepheid-like sources by inspecting the light curves with large when 
fitting a constant base-line. An example light curve is shown in Fig. [TTl These Cepheid-like 
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Fig. 7. — Same astrometric residuals as shown in Fig. [6] plotted as a function of position 
within the skycell. The gaps between the 4 detectors that overlap with this skycell are clearly 
visible (compare Fig. H]). 



- 22 - 



11.6 11.7 11.8 11.9 12.0 

RA [deg] 



41.0 41.1 41.2 

DEC [deg] 



l.Or 
C 0.5 

< 0.0- 

< -0.5 
-l.oL 



Fig. 8. — Absolute astrometric precision measured by comparing the positions of all sources 
in the rpi-band reference images of skycell 068 to a catalog of the Sloan Digital Sky Survey 
DR7. There is a small offset in Declination and a trend in Right Ascension as a function of 
Declination. Both effects are at the sub-pixel level (~ 0.1"). 



Skycell.077 PSF histrogiam Skycell.091 PSF histrogram 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 

PSF [arcsec] PSF [arcsec] 



Fig. 9. — Point-Spread Function (PSF) distribution of the first season of PAndromeda in the 
bulge (left panel) and disk (right panel) field. The PSF distributions of the rpi (ipi) data 
are shown in the blue (red). 



sources and further candidates that will be found in the remaining M31 fields will be analyzed 
in a forthcoming publication (Gossl et al. in prep.). In addition, we find one eclipsing binary 
and present its folded light curve in Fig. [T21 

During the time span of the first observing season of PAndromeda we find a few nova 
candidates in the central 40' x 40' sub-field. An example light curve is shown in Fig. [131 
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Fig. 10. — Percentage photometric accuracy (rms) as a function of the magnitude of stars 
from images taken by PAndromeda in the bulge field (left panel) and the disk field (right 
panel). The RMS is derived from the light curves of resolved stars. The black line is the 
expected total RMS, taking into account the contributions from sky (red line), read out 
(green line), object (blue line) and scintillation noise (magenta line). 
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Fig. 11. — Folded rpi-band light curve of our best Cepheid candidate at RA(J2000) = 
10.2662 deg and Dec(J2000) = 41.1512 deg. The period is 45.6 days and the distance from 
the center of M31 is 20'. 
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Fig. 12. — Folded rpi-band light curve of a foreground eclipsing binary at RA(J2000) 
10.6698 deg and Dec(J2000) = 41.1681 deg. The period is 0.5 day. 
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Fig. 13. — rpi-band light curve of the nova M31N-2010-10c from PAndromeda survey at 
RA(J2000) = 11.1108 deg and Dec(J2000) = 41.5205 deg. 
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5. The first microlensing events 

To compare with the previous resuhs, we first analyzed the central region of the M31 
field (40' X 40') and detect 6 candidates. The positions and light curves of these 6 candidates 
are shown in Fig. [TU Fig. [15] and Fig. [161 The fitting parameters (as shown in equation 
on])) and the x^of are presented in Table [S] 




Fig. 14. — Position of the six microlensing event candidates detected in the central 40' x 40' 
region of M31 from PAndromeda. The coordinates, RA (J2000) in hour and Dec (J2000) in 
degree are also shown in the figure. The FOV of this image is 40' x 40'. 



Table 3: PAndromeda microlensing event candidates in rpi - and zpi-band. 



lAU name 



RA(2000) 
[h:m:s] 



DEC(2000) 


^M31 


to 


^FWHM 


-^eff ,rpi 


-^cfE,ipi 


[d:m:s] 






[day] 


[ABmag] 


[ABmag] 


41:11:14 


10.0' 


5479.05±0.04 


1.28±0.09 


20.52±0.05 


20.29±0.06 


41:15:23 


4.2' 


5433.33±0.04 


2.07±0.16 


20.22±0.06 


19.80±0.08 


41:17:18 


2.3' 


5451.49±0.04 


3.34±0.14 


18.95±0.03 


17.60±0.05 


41:13:20 


9.5' 


5484.73±0.11 


13.82±0.43 


20.41±0.02 


20.07±0.03 


41:09:17 


19.3' 


5498.05±0.17 


9.04±0.36 


20.92±0.04 


20.18±0.05 


41:32:04 


27.8' 


5511.76±0.12 


3.08±0.40 


21.97±0.08 


21.27±0.12 



[ABmag] 



Xdof 



rpi 



PAnd-1 PSO JOlO.5746+41.1872 00:42:17.9 

PAnd-2 PSO JOlO.6329+41.2564 00:42:31.9 

PAnd-3 PSO JOlO.6596+41.2883 00:42:38.3 

PAnd-4 PSO JOlO.7979+41.2222 00:43:11.5 

PAnd-5 PSO JOlO.4579+41.1547 00:41:49.9 

PAnd-6 PSO JOlO.9700+41.5344 00:43:52.8 



0.23±0.10 
0.42±0.13 
1.35±0.08 



1.1 
1.2 
1.4 



0.34±0.06^ 1.1 
0.74±0.07^ 



0.69±0.22 



1.2 
1.0 



We give the lAU name of the 6 microlensing candidates in column 2, the coordinates in columns 3 and 4, the distance to the center of M31 
(Am3i) in column 5, the time of flux maximum (JD - 2450000) and event timescale in columns 6 and 7 the maximum flux in units of 
magnitude in rpi and ipi-band in columns 8 and 9, the color in column 10 and the best-fitted Xdof in columns 11 and 12. 



"The rpi and ipi are rather similar to the SDSS analogs fsee Tonrv et al. 2012 . for further description.). 



Table 4: PAndromeda microlensing event candidates with magnitudes in R and /-band. 









(R-I) 






[mag] 


[mag] 


[mag] 




PAnd-1 


20.10 


19.86 


0.24 




PAnd-2 


19.75 


19.32 


0.43 




PAnd-3 


18.21 


16.93 


1.28 


to 
-a 


PAnd-4 


19.96 


19.61 


0.35 




PAnd-5 


20.36 


19.63 


0.73 




PAnd-6 


21.42 


20.72 


0.70 





For readers who are used to the Johnson-Cousins system, we 
provide the magnutude and color in R and / from equation ([5]). 
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Fig. 15. — Light curves of the microlensing events detected in the central 40' x 40' region 
of M31 from PAndromeda. Upper (middle) panel shows the light curve in rpi (^pi) and the 
lower panel shows the residuals to the best-fitted light curve. In the figure there are also 
the name of the event (e.g. PAnd-1), the coordinates a (RA) and 5 (Dec) at the epoch 
of J2000, the distance to the center of M31 (Am3i) in arcminutes. The best-fit light curves 
(green) and parameters (black) are also shown, which are the time at maximum magnification 
(to), the event timescale (tpwHivi) in units of a day, the equivalent magnitude at maximum 
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Fig. 16. — Light curves of PAnd-5 and PAnd-6. See the figure caption of Fig. [15] for further 
explanation. 
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The 6 PAndromeda events form 3 groups. PAnd-2 and PAnd-3 are at 4.2' and 2.3' 
(0.5-1.0 kpc) from the Andromeda center. This region will be bulge-bulge lensing dominated 
regardless of the halo- MACHO content. The colors (rpi — zpi = 0.4 and 1.4 mag, R — I = 
0.4 and 1.3) and short time scales (tpwHM = 2.1 and 3.3 da y) of the events are in agreement 
with bulge stars lensing post main sequence bulge stars (see iRiffeser et a/.ll2006l ). The events 
PAnd-4 and PAnd-1 are at 9.5' and 10' distance from the M31 center along the minor and 
major axes (corresponding to ~ 2 kpc). These are regions where the light prof ile turns from 



bulge light dominated to disk light dominated (compare, e.g. the Fig. A.l of iTempel et al. 



201 ll . for the (?-band dust corrected minor and major M31 surface brightness profiles). Hence, 
the bulge-bulge and to a smaller degree the bulge-disk self-lensing rates are reduced relative 
to the center. On average the sources of self-lensing become more likely disk stars, and thus 
bluer. The timescales of self-lensing stay short, i.e. a few days. PAnd-1 has a timescale of 
1.3 days, the other (PAnd-4) is considerably longer and has a timescale of 14 days. This is 
ex pected only for a ve ry small fraction of self-lensing events (compare the density contours 
of 



Riffeser et a/.ll2006l . Fig. 9) and is more easy to reconcile with halo lensing. 



The events PAnd-5 and PAnd-6 finally are at 19.3' and 27.8' or 4.5 and 6 kpc from the 
M31 center along the major axis. This is the location where the SFB from the bulge drops 
to or below that of the young disk component, which is about 4 magnitudes below the SFB 
of the disk. At this distance from the M31 center bulge lensing should be almost zero and 
disk-disk self-lensing should be (because of the inefficient geometry due to the small line of 
sight extension) small as well. Halo lensing by MACHOs will decrease less fast out to 6 kpc 
since the disk stars provide a spatially almost 'flat' (compared to the bulge) resource for 
halo lensing at 2-6 kpc. The timescale of the two outer lensing events is 3 (PAnd-6) and 
9 (PAnd-5) days. The colors are both 0.7 mag (in either r-pi — i-pi or R — I), which is in 
agreement both with evolved bulge and disk stars. Since PAnd-5 is at a location where dust 
lanes become visible we also checked the extinction values along the line of sights to the 
events, to flnd out, whether one event could be pa rticularly reddened by M31 dust. We used 
the dust extinction map of iMontalto et all fl2009[ ) (the spatial resolution of the dust map is 
6"/pixel) and obtained E(B-V) values of 0.17 to 0.2 for all 6 events. According to this map 
PAnd-5 is close but not along the line of sight to a larger dust extinction. 

Our results concerning these long timescale events at 2 and 6 kpc from the Andromeda 
center are in agreeme nt with previous results: L ong time scale lensing events have been found 
by POINT- AGAPE dCalchi Novati et aZ-llioOsl ) at a distance of 22' from M31 (PA-99-N2) 
with a times scale of 22 days. 

We also inspect the HST images at the location of these events (see Fig. [T9|) . At this 
point we assume that the world coordinate system (WCS) of the HST images is correct and 
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we do not map them to the PSl images. There are HST archive images for all the events 
except PAnd-5. At the position of PAnd-3, PAnd-4 and PAnd-6, there are one to several 
resolved sources around the position of our candidates, while for PAnd-2 there is a possible 
source and for PAnd-1 nearby sources can be ruled out. 

We have not found extended sources (background galaxies) overlapping or close to our 
events, so we can rule out the contamination from SNe. This is also confirmed by the good 
Xdof from fitting the microlensing light curve. 
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Fig. 17. — Distribution of the -^cfr.rpi of the 6 PAndromeda microlensing events from Table 
|5l Most of the events have rpi fainter than 20 mag at the flux maximum except PAnd-3. 
Such a bright event is hard to reconcile with self-lensing due to the flux excess limit caused 
by finite-source effec ts. In general very bright events can be more easily caused by MACHOs 
faiffeser et oZlboOsI ). 
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Fig. 18. — Distribution of the tpwHM of the 6 PAndromeda microlensing events. The 
timescale of the detected events ranges from 1.3 to 13.8 days, with 4 out of the 6 events 
have tpwHM below 5 days. This demonstrates that the PAndromeda survey is effective in 
finding short duration events. 
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Fig. 19. — Comparison of PSl images (left column) with the HST images (right column). The 
HST image of PAnd-1 is from a median of 3 ACS 814 nm archive data (u3d9020at_drz.fits, 
u3d9020gt_drz.fits, u3d9020mt_drz.fits) of HST Cycle 5 proposal 6300 by Ford, H.. The HST 
image of PAnd-2, PAnd-3 and PAnd-4 are from ACS 814 nm archive data (jbf309010_drz.fits, 
jbf308010_drz.fits and jbf310010_drz.fits) of HST Cycle 18 proposal 12058 by Delcanton, J.. 
The HST image of PANd-6 is from ACS 606 nm archive data (j96q05010_drz.fits) of HST 
Cycle 13 proposal 10407 by Morrison, H.. 
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6. Outlook 

We are currently analyzing the full data set of the first PAndromeda season. Besides 
the point-source point-lens modeling and detecti on present e d in section [5l a search for the 



finite-source effects according to the process by iLee et al\ (120091 ) will also b e carried out. 



Furth er investigations on the individual bright events following the analysis of iRiffeser et al. 



( I2OO8I ) will shed light on the identity of these events, whether they are more likely caused 
by MACHOs than by stellar lenses. A comprehensive study of the detection efficiency (e.g. 
as a function of Ap and the event duration in the different fields) will also be done. 

In addition to the microlensing analysis, we will use the PAndromeda dat a-set to improve 



the M 31 model in several aspects, for example, including the dust content (IMontalto et al. 



20091 ). modeling the three-dimensional distribution of stars. At the time of writing this paper, 
we are also granted time to investigate the stellar population and dynamics of the M31 bulge 
in more detail (10 nights of the 2.7m telescope at McDonald observatory, with the VIRUS-W 
instrument), e.g. the bulge-bulge self-lensing (main contribution to the self-lensing in the 
center of M31) can be quantified at a much better level than now. 

The IPP is going to deliver a reprocessing of the PAndromeda data in some months 
from now. We will analyze the reprocessed data and also the upcoming new data without 
binning. The observation strategy in 2011 is similar to 2010. We have two rotation angles of 
the de-rotator both with a good matching on a grid of corresponding detectors. In 2010 the 
stars on the central part of the cameras were defocused, leading to a decreased S/N. This 
can mimic an asymmetric microlensing signal. Therefore for 2011 we shift this defocused 
center to the major axis of M31 to have more symmetric signals from both sides of M31 (see 
Fig. [T]). Of course, our final detection efficiency study will have to account for the varying 
PSF sizes in different skycells. 

In 2011 we continue with our two block strategy with a time gap of three to five hours 
accepting short offsets of <30 min due to scheduling constraints. The filters and exposure 
times will be the same like in the previous season. 

The pointing follows a dither pattern with a diameter of 13" to 20" and a cycle of 5 
images. This ensures that all the detector gaps are filled completely. The regular observations 
have started on 25th July 2011. 
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7. Conclusion 

The preliminary analysis of the first part of the first season PAndromeda data demon- 
strates that one can indeed detect microlensing events at a competitive rate per season (see 



Calchi Novatil |2010| . for a summary of previous surveys and events) . The data are further 



useful for novae detection and investigation of other variables. 

The identification of 4 short-duration microlensing events with tpwHM ~ 1-3 days shows 
that the time resolution of the PAndromeda project is comparable with the best 2 seasons 
of the WeCAPP project (where two telescopes were coordinated to monitor M31). The 
advantage of PAndromeda is the capability to cover the entire M31 disk with one pointing 
(compared to the 17.2' x 17.2' FOV of WeCAPP). By finding or excluding microlensing 
events in the outer disk of M31, the PAndromeda survey will manifestly exceed the accuracy 
of previously derived M31-halo MACHO fractions. 
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